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Abstract 

The method for calculation of valence quark distributions at intermediate x is pre- 
sented. The imaginary part of the virtual photon forward scattering amplitude on the 
quark current with meson quantum number is considered. Initial and final virtualities 
p\ and p\ of the currents are assumed to be large, negative and different, p\ ^ p\. 
The operator product expansion (OPE) in p\,p\ U P to dimension 6 operators is per- 
formed. Double dispersion representations mpi,p| of the amplitude in terms of physical 
states contributions are used. Equalling them to those calculated in QCD by OPE the 
desired sum rules for quark distributions in mesons are found. The double Borel trans- 
formations are applied to the sum rules, killing non-diagonal transition terms, which 
deteriorated the accuracy in the previous calculations of quark distributions in nucleon. 
Leading order perturbative corrections are accounted. Valence quark distributions in 
pion, longitudinally and transversally polarized p-mesons are calculated at intermediate 
x, 0.2 < x < 0.7 and normalization points Q 2 = 2 — 4 GeV 2 with no fitting parameters. 
The use of the Regge behaviour at small x and quark counting rules at large x allows 
one to find the first and the second moments of valence quark distributions. The ob- 
tained quark distributions may be used as an input for evolution equations. In the case 
of pion the quark distribution is in agreement with those found from the data on the 
Drell-Yan process. Quark distributions in transversally and longitudinally polarized 
p-mesons are essentially different. 
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1 Introduction 



Quark and gluon distributions in hadrons are not calculated in QCD starting from the first 
principles. What is only possible is to calculate their evolution with Q 2 . In the case of 
nucleon the standard way is the following. At some fixed Q 2 = Ql the form of quark and 
gluon distributions, characterized by the number of free parameters is assumed. Then the 
evolution of distributions with Q 2 is calculated and after comparing with the data of the deep 
inelastic lepton-nucleon scattering the best fit for these parameters is found. (About dozen 
of the parameters is used, see, e.g. the recent papers [0-01). Quark and gluon distributions 
in pions are determined in a similar way from the data on production of e + e~ and pairs 
in pion- nucleon collisions (the Drell-Yan process) || § and prompt photons production (?]]. 
But here one needs an additional hypothesis about connection of the fragmentation function 
at time-like Q 2 , which is measured in the Drell-Yan processes, with parton distributions 
defined at space-like Q 2 . For other mesons and baryons, as well as for unmeasured up to 
now distributions in polarized nucleon (like h\{x)) we have no information about parton 
distributions from experiment and have to rely on models. The QCD calculation of quark 
and gluon distributions in hadrons, which could be used as an input in evolution equations 
may be considered challenge for the theory. 

The method of valence quark distributions calculation at intermediate x was suggested in 
H and developed in p|-pT||. The idea was to consider the imaginary part (in s-channel) of 4- 
point correlator, corresponding to the forward scattering of two-quark currents, one of which 
has the quantum numbers of hadron of interest and the other is electromagnetic (or weak). 
It was supposed that virtualities of the photon and hadronic current q 2 and p 2 are large and 
negative \q 2 \ 3> \p 2 \ 3> R^ 2 , where R c is the confinement radius (q is the momentum of 
virtual photon, p is the momentum, carried by hadronic current). It was shown ||, that in 
this case the imaginary part in s channel [s = {p + q) 2 } of the forward scattering amplitude is 
dominated by small distance contributions at intermediate x. (The standard notation is used: 
x is the Bjorken scaling variable , x = —q 2 /2i>, v = pq). The proof of this statement, given 
in H, is based on the fact that for the imaginary part of the forward scattering amplitude 
the position of the closest to zero singularity in momentum transfer is determined by the 
boundary of the Mandelstam spectral function and is given by the equation 



Therefore, if \p 2 \ is large and x is not small, then even at t = (the forward amplitude) the 
virtualities of intermediate states in t-channel are large enough for OPE to be applicable. The 
further procedure is common for the QCD sum rules. On one side, the amplitude is calculated 
by OPE with the account of condensates. On the other side, the dispersion representation 
in p 2 in terms of physical states is written for the same amplitude and the contribution of 
the lowest state is extracted by using the Borel transformation. As follows from (HI), the 
approach is invalid at small x. (Since in real calculations — p 2 are of order 1 — 2 GeV 2 , small 
x means in fact x < 0.1 — 0.2). This statement is evident apriori, because at small x Regge 
behaviour is expected, which cannot be described in the framework of OPE. The approach 
is also invalid at x close to 1. This is the domain of resonances, also outside the scope of 
OPE. The fact that this method of calculation of quark distributions in hadrons is invalid 
at i C 1 and at 1 — x <C 1 follows from the theory itself: the OPE diverges in these two 
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domains. Therefore, calculating higher order terms of OPE makes it possible to estimate up 
to which numerical values of x, in the small and large x domain, the theory is reliable in each 
particular case. In the way described above, valence quark distributions in nucleon || |TTJ 
were calculated. However, the accuracy of the calculation was not good enough, especially 
for (i-quarks ||. Moreover, it was found to be impossible to calculate quark distributions in 
7r- and p-mesons in this way. The reason is that the sum rules in the form used in [TTJ 



have a serious drawback. (The calculation of the photon structure function |T(J is a special 
case and has no such problem). 

The origin of this drawback comes from the fact that in the case of 4-point function, which 
corresponds to the forward scattering amplitude with equal inital and final hadron momenta, 
the Borel transformation does not provide suppression of all excited state contributions: the 
non-diagonal matrix elements, like 

(0\j h \h*)(h*\f l (x)f l (0)h)(h\j h \0) (2) 
are not suppressed in comparison with the matrix element of interest 

(Ol/l/iX/il/^COJI/iX/il/lO) (3) 

proportional to the hadron h structure function. (Here h is the hadron which structure func- 
tion we would like to calculate, h* is the excited state with the same quantum numbers as 
h, j h is the quark current with quantum numbers of hadron h, j el is the electromagnetic 
current). In order to kill the background matrix elements @ it was necessary to differentiate 
the sum rule over the Borel parameter. But, as is well known, differentiation of an approxi- 
mate relation may seriously deteriorate the accuracy of the results. In QCD sum rules such 
differentiation increases contributions of higher order terms of OPE and excited states in 
physical spectrum, the sum rules become much worse or even fail (as for 7r and p- mesons). 
For pion the situation is especially bad, because direct calculation show, that the leading 
term in OPE (the bare loop diagram) corresponds just to the non-diagonal matrix element, 
not to the pion structure function. 

In ref.12 it was suggested the modified method of calculation of the hadron structure func- 
tions (quark distributions in hadrons), where this problem is eliminated and valence quark 
distributions in pion were calculated. This method is used here for calculation of valence 
quark distributions in p-meson, separately for longitudinally and transversally (relative to 
the virtual photon beam) polarized p-mesons. Since quark distributions in p-meson cannot 
be measured experimentally, the common way to get them is to assume SU(6) symmetry, 
where ir- and p-mesons belong to the same multiplet. Then quark distributions in p-meson 
are equal to that in pion and are independent on p polarization. On the other side, the pion 
play a specific role in the theory - it is a Goldstone boson. From this point of view it has 
nothing in common with p-meson and there are no reasons to expect, that quark distributions 
in 7r and p are equal. This problem will be resolved and it will be shown, that valence quark 
distributions in pion, longitudinally and transversally polarized p-mesons are quite different. 



For reader convenience we present first shortly the method of [T2| and the results for 
quark distributions in pion. Then valence quark distributions in longitudinal and transverse 
p mesons are calculated. The valence quark distributions are reliably calculated in the domain 
of intermediate x, 0.2 < x < 0.65. The use of the Regge behaviour at small x and quark 
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counting rules at large x allows one to find the first and the second moments of quark 
distributions. 



2 The method 

Consider the non- forward 4-point correlator: 

n(pi,p 2 ;g,g') = - 



d xd yd ze 



ip\x+iqy—ip2Z 



x (0|T{/(x), f(y), f(0), j h (z)}\0) 



(4) 



Here p\ and P2 are the initial and final momenta carried by hadronic current j , q and 
q' = q + p±—p2 are the initial and final momenta carried by virtual photons. (Lorenz indeces 
are omitted). It will be very essential for us to consider non-equal pi,p2 and treat p\,p\ as 
two independent variables. However, we may put q 2 = q' 2 = q 2 and t — (pi — P2) 2 = 0. We 
are interested in imaginary part of T\.(p\,p\, q 2 , s) in s channel: 



ImU(pl,pl,q 2 ,s) = — 



R(p 2 i,pl,Q 2 ,s + ie) -Yl(p 2 ,pl,q 2 ,s-ie) 



(5) 



In order to construct representation of ImU(p 2 ,p2, q 2 , s) in terms of contributions of physical 
states, let us write for Imll^^p 2 ,, q 2 , s) the double dispersion relation in p\,p\: 

¥(q 2 ,s,u) 7p(q 2 ,s,u) 



ImYl(plplq 2 ,s) = a(q 2 ,s) + f ™ ' % Mj du + f f^L 

J u — pi J u — 



-du 



+ / dui / du 



p(q 2 ,s,u 1 ,u 2 ) 
(u x -pj)(u 2 -pl) 



(6) 



The second and the third terms in the right-hand side (rhs) of (§) may be considered as 
subtraction terms to the last one - the properly double spectral representation. The first 
term in the rhs of (§) is the subtraction term to the second and third ones. Therefore, (f|) has 
the general form of the double spectral representation with one subtraction in both variables 
- p\ and p\. Apply the double Borel transformation in p\,p\ to @. This transformation kills 
three first terms in rhs of (161) and we have 



® Ml® Ml Iml\.(p\,p 2 2 ,q 2 ,s) 



du\ / du2p{q , s, Mi, U2)exp 



Ux 



U2_ 

M 2 



(7) 



The integration region over U\,u 2 may be divided into 4 areas (Fig.l): 

I. ui < s ; u 2 < s 

II. Mi < s ; u 2 > so 

III. Mi > s ; m 2 < s 

IV. Mi,M 2 > s 

Using the standard QCD sum rule model of hadronic spectrum and the hypothesis of 
quark-hadron duality, i.e. the model with one lowest resonance plus continuum, one may 
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clearly see, that area I corresponds to resonance contribution. Spectral density in this area 
can be written as 

p(u 1 ,u 2 ,x,Q 2 ) = g 2 h ■ 2nF 2 (x,Q 2 )5(u 1 - m 2 h )5{u 2 - m 2 h ), (8) 
where is defined as 

(0\jh\h)=g h (9) 

(For simplicity we consider the case of the Lorenz scalar hadronic current. The necessary 
modifications for cases of tt and p-mesons will be presented below). If in ImU(pi,p 2 ,q,q') 
the structure, proportional to P^P U [P M = (pi +p 2 ), Jl /2] is considered, then in the lowest twist 
approximation F 2 (x, Q 2 ) is the structure function depending on the Bjorken scaling variable 
x and weakly on Q 2 = —q 2 . 

In area (IV), where both variables U\ y2 are far from resonance region, the non-perturbative 
effects may be neglected, and as usual in sum rules, the spectral function of hadron state is 
described by the bare loop spectral function p° in the same region 

p(ui,u 2 ,x) = p°(u 1 ,u 2 ,x) (10) 

In areas (II), (III) one of the variables is far from the resonance region, but other is in the 
resonance region, and the spectral function in this region is some unknown function p = 
ip(ui,u 2 ,x), which corresponds to transitions like h — ► continuum as shown in Fig. 2. After 
double Borel transformation the physical side of the sum rule can be written as (M 2 , Mf are 
Borel mass square) 



B x B 2 [ImIl] = 2nF 2 (x, Q 2 ) ■ g\e 



i 

Mi 



+ / du\ 



du 2 ip(ui, u 2 , x)e M i 



*•<> 



+ / du x / du 2 i/;(ux,u 2 ,x)e ^ ^ + / du 1 du 2 p°{u ll u 2) x)e M V (11) 

so so 



In what follows we put M\ = Mf = 2M 2 . (As was shown in |T3|, the values of Borel 
parameters Mf , Mf in the double Borel transformation are about twice of that in the ordinary 
ones) . 

One of advantages of this method is that after double Borel transformation unknown 
contribution of (II), (III) areas [the second and third term in (|TT1) 1 are exponentially sup- 
pressed. Using duality arguments, we estimate the contribution of all non-resonancee region 
(i.e. areas II, III, IV) as a contribution of bare loop in the same region and demand their 
value to be small (less than 30%). So, equating physical and QCD representation of II and 
taking into account cancellation of appropriate parts in the left and right sides, one can write 
the following sum rules: 

_ m 2 ( 1 j 1 ) 

Im U° QCD + Power correction = 2ttF 2 {x, Q 2 )g 2 h e h 

SO SO 

ImU° QCD = J J p°(ux,u 2 ,x)e-^ (12) 
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It can be shown (see below), that for box diagram i/;(ui,U2,x)) ~ 5{u\ — U2), and, as a 
consequence, the second and third terms in ( |Tl~D are zero in our model of hadronic spectrum. 

It is worth mentioning that if we would consider the forward scattering amplitude from 
the beginning, put Pi = P2 = P and perform Borel transformation in p 2 , then unlike (]TT|), the 
contributions of the second and third terms in (Q) would not be suppressed comparing with 
the interesting for us lowest resonance contribution. They just correspond to the non-diagonal 
transition matrix elements discussed in the Introduction and are proportional to 



(0\j h \h*) -^—^ (h*\f l (x)f\0)\h) -jJ-^fcl/lO) (13) 
p z — mjf p z — m^ 



From decomposition 



p 2 — m* 2 p 2 — m\ m* 2 — m 2 h \p 2 — m* 2 p 2 — m 2 h 



(14) 



it is clear that in this case (|13|) may contribute to the second (or third) term in (|6|) and 
after Borel transformation the contribution of the second term in flT4] ) has the same Borel 
exponent e~ m ^l M as the lowest resonance contribution. The only difference is in pre-exponent 
factors: they are 1/M 2 in front of the resonance term and Const, in front of non-diagonal 
terms. This difference was used in order to get rid of non-diagonal terms: application of 
the differential operator (d/d(l/M 2 )e m h/ M to the sum rule kills the Borel non-suppressed 
nondiagonal terms, but deteriorates the accuracy and shrinks the applicability domain of the 
sum rule (particularly, the domain in x, where the sum rule is valid). 

Show now that in the used here model of hadronic spectrum-resonance plus continuum, 
where continuum is given by the bare loop contribution, the second and third terms in (|TT|) 
are in fact zero. Consider the bare loop represented by the diagram of Fig. 3. For simplicity 
restrict ourselves by the case when all propagators are bosonic and all currents are scalar. (In 
the realistic case with quarks in internal lines conclusion will be the same). The imaginary 
part in s-channel of Fig.3 diagram is given by (quark masses are neglected): 



Im T(plp 2 2 , q 2 ,s) = J *k {pi _ k) t {p2 _ k)2 6 



(Pi + Ql ~ kf 



8{k 2 ) (15) 



Neglecting the higher twist terms ~ Pi 2 /<? 2 (HI) is equal to (see [[TI|, Appendix) 

Im T(p 2 v pj, q 2 , s) = - — / 2-- 2T du ( 16 ) 

4z/x J [u — Pi){U — P2) 



where v = qP and x = —q /2u. To derive (It) it is convenient to introduce 



P= (Pi+p 2 )/2 

r — Pi — Vii r 2 = (17) 

and in the Lorenz system, where 4-vector P has only ^-component and q - only time and 
z-components 

r °~ rz ~2^' r± ~ 4 v 2 -q 2 P 2 ' 
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r\ is of higher twist and may be neglected. As follows from flTB|), in case of bare loop of 
Fig.3 the spectral function is proportional to 8(ui — w 2 ) and the contributions of areas II, III 
in Fig.l are zero. Since one may expect, that more complicated diagrams are smaller and in 
any case their contributions in domains II, III are suppressed by Borel exponents, one may 
safely neglect them, as was done in (|T^). 



3 Quark distributions in pion 

It is enough to find the distribution of valence w-quark in 7r + , since d(x) = u(x). The most 
suitable hadronic current in this case is the axial current 



j»5 = u'y^d (19) 

In order to find the w-quark distribution, the electromagnetic current is chosen as w-quark 
current with the unit charge 



J, = ui,u (20) 
The bare loop Fig.3 contribution is given by 

Im U ^ = ~ \ I if (k+pl-p^ + k)2]5[pi ~ k)2] 

x Tr \j\k-f ll (k + q)^ u (k +p2 - Pi)la(k - pi)} (21) 

The tensor structure, chosen to construct the sum rule is a structure proportional to 
P^PyPxPfj/v. The reasons are the following. As is known, the results of the QCD sum 
rules calculations are more reliable, if invariant amplitude at kinematical structure with 
maximal dimension is used. Different pi ^ P2 are important for us only in denominators, 
where they allow one to separate the terms in dispersion relations. In numerators one may 
restrict oneself to consideration of terms proportional to 4-vector P^ and ignore the terms 
~ r M . Then the factor P^Py provides the contribution of i^x) structure function and the 
factor P\P a corresponds to contribution of spin zero states. (The factor \ jv is scaling 
factor : W2 = F 2 /v.) 

Let us use the notation 

IW = {P»PvPxPa/v)Il(p 2 i,pi x) + ... (22) 

Then ImU(pl,p2, x) can be calculated from (|21|) (eq.'s (|T5|),(|T6|) are exploited) and the result 
isH: 

Imll(pl,pl,x) = -x 2 (l-x) [ dui [ du 2 -, — ^ — ^ (23) 
Tr oo (^i -Pi)( M 2 -P2) 

The matrix element of the axial current between vacuum and pion state is well known 

(0\j^\n)=tpJ 7T (24) 

where f n = 131MeV. The use of (0), (|23D , (|24]) gives the sum rule for valence w-quark 
distribution in pion in the bare loop approximation |T2|j : 
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= A " x ^ ~ e-^ M2 )e</ M \ (25) 

J TT 



where Sq is the continuum threshold. In ref. |L2[ the following corrections to (P5|) were ac- 
counted: 

1. Leading order (LO) perturbative corrections, proportional to ln(Q 2 /fi 2 ) , where yii 2 is 
the normalization point. In what follows the normalization point will be chosen to be equal 
to the Borel parameter /i 2 = M 2 . 

2. Power corrections - higher order terms of OPE. Among the latter, the dimension- 
4 correction, proportional to gluon condensate {0\^-G™ u G™\0) was first accounted, but it 
was found that the gluon condensate contribution to the sum rule vanishes after double 
borelization. There are two types of vacuum expectation values (v.e.v) of dimension 6: one, 
where only gluonic fields enter: 

^a s f abc (0\G; u G b vX Gyo) (26) 
and the other, proportional to four-quark operators 

<o|^r-0 - -0r-0|o> (27) 

It was shown in ]T^] that terms of the first type cancel in the sum rule and only terms of the 
second type survive. For the latter one may use the factorization hypothesis which reduces 
all the terms of this type to the square of quark condensate. 

A remark is in order here. As was mentioneed in the Introduction, the present approach 
is invalid at small and large x. No-loop 4-quark condensate contributions, like Fig.4, are 
proportional to 8(1 — x) and being outside of the applicability domain of the approach, 
cannot be accounted. In the same way, the diagrams, which can be considered as a radiative 
corrections to those, proportional to 8(1 — x), must be also omitted. 

All dimension-6 power corrections to the sum rule were calculated in ref. 12 and the final 
result is given by (the pion mass is neglected): 

J^M 2 2 Jf^ , (a s (M 2 )-ln(Q 2 /M 2 



X 



1 + Axln(l — x) 




x 

AiTa s (M 2 ) ■ 4ira s a 2 uj(x) 



(2tt) 4 • 3 7 • 2 6 • M 6 x 3 (l-x) 3 
where uj(x) is 4-order polynomial in x, 



(2? 



a = -(2tt) 2 (0|^|0) (29) 

uj(x) = -5784x 4 - 1140x 3 - 20196a; 2 
+20628x - 8292)/n(2) + 4740x 4 + 8847x 3 

+ 2066x 2 - 2553x + 1416 (30) 
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u n (x) may be used as an initial condition at Q 2 = Qq for solution of QCD evolution equations 
(DGLAP) equations Q. 

In numerical calculations we choose: the effective LO QCD parameter Ag^. D = 
200AfeV, Ql = 2GeV 2 } a s a 2 (lGeV 2 ) = 0.13GeV 6 |[J. The continuum threshold was 
varied in the interval 0.8 < So < 1.2GeV 2 and it was found, that the results only slightly 
depend on it. The analysis of the sum rule ( P5| ) shows, that it is fulfilled in the region 
0.15 < x < 0.7; the power corrections are less than 30% and the continuum contribution 
is small (< 25%). The stability in the Borel mass parameter M 2 dependence in the region 
OAGeV 2 < M 2 < 0.6GeV 2 is good. The result of our calculation of valence distribution in 
pion xu n (x, Ql) is shown in Fig. 5. 

Suppose, that at small x < 0.15 u n (x) ~ 1/v^ according to Regge behaviour and at 
large x > 0.7 u n (x) ~ (1 — x) 2 according to quark counting rules. Then, matching these 
functions with (p8|), one may find the numerical values of the first and second moments of 
w-quark distribution 

i 

Mi = J u n (x)dx « 0.84 (0.85) (31) 
o 
i 

M 2 = J xu w (x)dx ss 0.21 (0.23) (32) 
o 

(In the parenthesis we give the values, corresponding to u w (x) ~ 1 — x behaviour at large 
x.) The results only slightly depend on the points of matching (not more than 5%, when 
the lower matching point is varied in the region 0.15 - 0.2 and the upper one in the region 
0.65 - 0.75). M.\ has the meaning of the number of w-quarks in ir + and should be Ai = 1. 
The deviation of ([H]) from 1 characterizes the accuracy of our calculation. Ai 2 has the 
meaning of the part of pion momentum carried by valence w-quark. Therefore, valence u and 
d quarks are carrying about 40% of the total momentum. In Fig. 5 we plot also the valence 
it-quark distribution found in || by fitting the data on production of yU + yU~ and e + e~ pairs 
in pion-nucleon collisions (Drell-Yan process). Comparing with the found here distribution 
it must be taken in mind, that the accuracy of our curve is of order of 10 — 20%, the last 
number refers to the border of the applicability domain. ?7-quark distribution found from 
the data on the Drell-Yan process is also not free from uncertainties. Particularly, what 
is measured in the Drell-Yan process is the quark fragmentation function into pion defined 
at q 2 > 0. In order to get quark distribution in pion defined at q 2 < 0, the procedure of 
analytical continuation is used, which may introduce some uncertainties, especially, at low 
normalization point, like Ql = 2 GeV 2 , to which the data in Fig. 5 refer. For all these 
reasons we consider the agreement of two curves as good. The calculation of valence w-quark 



distribution in pion in the instanton model was done recently fl4 |. At intermediate x the 



values of xu n {x) found in are about 20% higher, in comparison with ours. Recently also 
the pions valence quark momentum distribution using a Dyson-Schwinger equation model 
was found in JT5fl . Our results are in reasonable agreement with the results of [15]. 



1 There was a misprint in the corresponding equation in |l2| (eq.[t(] of |[2]]): instead of a s (M 2 ) in the last 
term was a s (Ql). In numerical calculations the correct value a s (hP) was taken. 



9 



4 Quark distributions in p-meson 



Let us calculate valence w-quarks distribution in p + -meson (p - width is neglected). The 
choice of hadronic current is evident 

jfl = u ltl d (33) 
The matrix element (p + \ j p p | 0) is given by 

2 

Tfi 

(P + I % I 0) = (34) 
9p 

where m p is the p-meson mass, g p is the p — 7 coupling constant, g p /4:ii = 1.27, is 
the p-meson polarization vector. Consider the coordinate system, where the collision of p- 
meson with momentum p and virtual photon with momentum q proceeds along z-axes. The 
averaging over p polarizations is given by the formulae: 
1. Longitudinally polarized p: 

y m p J \ m p J v A — q A m p 



2. Transversally polarized p: 



^ ^ \ p m 2 p J u 2 - q 2 m 2 \" mj J \ mj J 

The imaginary part of the forward p — 7 scattering amplitude W^ u x a (before multiplication 
by p-polarizations) satisfies the equations: W^xaQ^ = W^xaQv = W^ vXa p x = W^XaPa = 0, 
which follows from current conservation. (The indeces p, v refer to initial and final photon; 
A, a - to initial and final p.) The general form of W^xa is: 



ft - If) [» - -#) {>>■ - ^rf + ^ - - -7- j 



« f ^ n 



(37) 

where A,B,C,D are invariant functions. Average (|37D over polarizations for longitudinal and 
transverse p-mesons. We have 



/ 2 q 1 - 
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~ E W^e\el = - L v - q -f-y+(p»-^)[ P »-^)C (39) 




From comparison of ([!7]) and (|38|), ( j39|) it is clear, that the proportional to Pp,p v structure 
function Fi(x) in the scaling limit [y 2 q 2 \ m 2 ) is given by the contribution of invariants 
C + (u 2 /m 2 p )D in case of longitudinal and by invariant C in case of transversal p-mesons. This 
means, that in the forward scattering amplitude W^ v \ a ([37]) one must separate the structure 
proportional to p^p u P\Pa in the first case and the structure ~ p^Pv&Xa m the second case. 
Consider now the non-forward 4-point correlator 



^uXa(Pi,P2]q,q) = -if d A xd A yd A ze 



ipix+iqy—ip2Z 



x (0 | T {f x (x), #(y), fj(0), #(*)} | 0), (40) 

where the currents j^{x) and j\(x) are given by (|20|) and ([33|). It is evident from the 
said above, that in the non-forward amplitude for determination of w-quark distribution in 
longitudinal p-meson the most suitable tensor structure is that, proportional to P^PuPaPx, 
while M-quark distribution in transverse p can be found by considering the invariant function 
at the structure —P^P u S\ a (P is given by fllTf) ). 

Calculate first the bare loop contribution (diagram of Fig. 3). In case of longitudinal p- 
meson the tensor structure, which is separated is the same, as in the case of pion. Since 
at m q = bare loop contributions for vector and axial hadronic currents are equal, the 
only difference from the pion case is in the normalization. It can be shown, that w-quark 
distribution in longitudinal p-meson can be found from ( |25| ) by substitutions — > m p , 
/tt — > m p /g p and therefore 



3 M 

The calculation of g 2 performed in the same approximation in [To] leads to 



L p = 7^—9 2 Al - x)e^' M \l - e- s °' M2 ) (41) 



,2 71^-2 



m 2 4ir 

(the definition of g p used here differs from that in |l6j by a factor 1/V2). The substitution 
of © in (H) gives 

u L p (x) = 6x(l - x) (43) 

and we have 

l 

u L p (x)dx = 1, (44) 



o 



as it should be. Also, 



i 



o 



xu p (x)dx = -, (45) 
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is in correspondence with naive quark model, where w-quarks are carrying one half of to- 
tal momentum. The calculation of the term, proportional to the structure P^P v 5\ a in the 
diagram of Fig. 3 gives 



3 1 
2ttu 



— x(l — x) 



udu 



(u - pl)(u - p|) 



(46) 



After borelization we get for w-quark distribution in transversally polarized p-meson in bare 
loop approximation 



x(l — x) 

2 v ; 



(47) 



where 



E 1 {z) = 1- (l + z)e~ 2 (48) 

Let us use fl42|), put M 2 = m 2 and neglect the terms ~ e - s o/ m p_ Then a simple formula for 
v?(x) follows: 



x(l — x) 

2 v ; 



u p (x) = 3 

CZ-quark distribution ( f49l) has the expected properties: 



(49) 



u 1 p{x)dx = 1 



xu 1 p{x)dx 



(50) 



(51) 



Take in account LO perturbative correction, proportinal to In Q 2 j '/i 2 and choose Q 2 = Ql 
for the point where we calculate our sum rules. The result is (the second term in square 
brackets corresponds to the perturbative correction): 



and 



L,., - 3M2 9 L m V M2 x(l-x) 



<{x) 



47r 2 m 2 



1 a s (fx 2 )ln(Q 2 //i 2 ^ 
2(1 - 2x)ln x 



x [l/x + 4/n(l - x) - 



1 — x 



3tt 



(52) 



Sir \ 



2(1 - 2x + Ax 2 )lnx s 
(p (x) 



(53) 
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where 



(p (x) = 1 - 2x(l - x) (54) 

Turn now to consideration of power correction contribution to the sum rules. The power 
correction of lower dimension is proportional to the gluon condensate (G^G^) with d = 4. 
As was discussed above, only s-shannel diagrams (Fig.l) exist in the case of double bore- 
lization. The (G^G^) correction was calculated in a standard way in the Fock-Schwinger 
gauge x^Ap = pj]. 

The quark propagator iS(x, y) = (ip(x)tp(y)) in the external field Ap has the well-known 
form [13], H |i~8f . In the 7r-meson case the sum of all diagrams, corresponding to (G a G a ) 
corrections (Fig.6), was found to be zero after double borelization |12| . For p-meson, however, 



(G^G'L,) correction are equal to zero for longitudinally polarized p (pl) but are not vanishing 
for transversally polarized p(pr) 

Mr 4) = ~(0 I ^G\ v I 0) (55) 

All diagrams here and in what follows are calculated using a program for analytical calcu- 
lations REDUCE. Before we discuss d = 6 contribution, let us remind, that we should not 
take into account non-loop diagrams and diagrams, which can be treated as their evolution 
(for example see Fig. 7). There is a large number of loop diagrams for d = 6 correction. It is 
convenient to divide them into two types and discuss these types separately: 
Type I-diagrams where only interaction with the external gluon field is taken into account - 
see Fig. 8, 9. 

Type II - diagrams, where expansion of quark field is also accounted (V-covariant derivate) 

i>(x) = if>(0) + x ai [V Q1 V(0)] + ^x ai x a2 [V Q1 V a2 ^(0)] + ... (56) 

The examples of such diagrams are shown in Fig.10,11. Discuss briefly the special features 
of calculation of this two types of diagrams. The diagram of type I (Fig.8) are, obviously, 
proportional to (0 | g 3 j^G a G\^G C \ 0) and when calculating it is convenient to use the 



representation of this tensor structure suggested in [[19 



(0 | g 3 f abc G%G b a pG c pa | 0) = 1/24(0 | g 3 f abc G a yS G b Se G c ei \ 0) (57) 

X {giio9a.v9flp + 9iifi9a P 9cjv + 9ao9{ip9vP + 9pu9pa9p<r 



-gp,p9ao9 P u — g^gupgap — gavg^pgpa — g/3 P gp, a g, 



The diagrams of Fig.9 are proportional to (0 | D p Gl u D T G a ap \ 0) and (0 | G a pu D p D T G a a/3 | 0). 
Using the equation of motion it was found in Jl9| thatQ 



■(0 | D p G a pp D a G a a p | 0) = (0 | G a pv D P D a G a aS3 \ 0) = 2Q~ [g^g^ - g m g u p) 



2 Note that our choice of sign of g is opposite to those of [pl 
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J T'^\9p,fi9a.a9pv + 9uv9pp9po — 9aa9pp9vfi ~ 9pv9pa9fia)\ 

uv gpagapgvfi g pflg p,a.gi>(r) i 

O ± = ^(0 I 9 2 JX I °) T ^(0 I 9f abc G%G b vX G c Xll | 0), (58) 

where j« = E^7/,(A72)^. 

From (57) and (|58|) one may note that these tensor structures are proportional to two 



vacuum averages: 

(0\g%\0) and (0 | g 3 G a ^G b up G c pfl f abc | 0). 

The first of these, (0 | g 2 j 2 | 0), by use of the factorization hypothesis easily reduces to 
(gipip) 2 which is well known, 

(0|^|0) = -(4/3)[(0|^|0)] 2 . (59) 

But (0 | g 3 G pu G b /p G pfl f abc | 0) is not well known; there are only some estimates based on the 
instanton model [|20| ^TJ. In the 7r-meson case the terms, proportional to 



(0 | g z f abc G a pv G h vp G° | 0) exactly cancelled fll2"| . The similar cancellation takes place for pi. 



But there is no such cancellation for px and one should estimate (0 | g f G a ^ u G up G c pp | 0). 
The estimation, based on the instanton model [GO], gives 



- (9 3 f abc G%Gl p G c pp ) = -=--o<0 | (a s /n)Gl | 0), (60) 



487T 2 1 

where p c is the effective instanton radius. 

Among the diagrams of type II (Figs. 10, 11) only those, where the interaction with the 
vacuum takes place inside the loop, are considered. Such diagrams cannot be treated as the 
evolution of any non-loop diagrams and are pure power corrections of dimension 6. All these 
diagrams are, obviously, proportional to 

(0\^ b D p G; u \0), (0|^(V T $)G£,|0), (0 i (vw^fei 0). 

These tensor structure were considered in where with the help of equation of motion the 
following results were obtained^: 



9(0 | H> | 0) 



(0 | rMD a G, v ) n | 0) = 3 3 . 2 5 (fln/T/x - 9.plu)p a (X n r\ (61) 



(0 | ^ d a (V^) b G; u | 0) = 9( °| 3 # 2 | ° y [g^u - g*vj» - ^A757AK(A"r. (62) 

The term (0 | i^J a ^)a) d ^G n pv | 0) can easily be calculated using (|B~1"D , fl62"|) 
3 The sign errors in front of g in eq.'s (p|),(|3),(Al),(A2) of [111 are corrected. 
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(0 I (V^^G* I 0) 



9(0 I W I 0) ; 
3 3 -2 6 



For diagrams in Fig. 11 we use the following expansion of the gluon propagator: 



(63) 



S%{x-y,y) 



(27T) 



zip/ 



-ik(x-y) 



-^(kH aX -2k a k x ))(D a G px ) 



-ik x y a G l aX - -i(y a ypk x 



1 / 2iA; c 



<^p + 2 



(64) 



This expression can be found by the method of the calculation of the gluon propagator in 



the external vacuum gluon field suggested in |T7| (see also |TJ, |22|). The total number of 
d = 6 diagrams is enormous - about 500. Collecting the results together we get finally the 



following sum rules for u p 



XUp(x) 



An 2 



M 2 ^-e m "' M2 x 2 (l-x) 
mi 



1 + 



'a s {M 2 )-ln{Q 2 /M 2 
3vr 



x 



'l + Axln(l -x) 2(l-2x)lnx" 



x 



1 — x 



'l-e 



a s {M 2 



CUM 



Uj(x) 



71 



2 -3 7 -2 6 -M 6 x 3 (l-x) 



(65) 



where a and u(x) are given by ( p9|) , fl3~0"|). A remark here is in order. In (|65|) (as well in 
the strong coupling constant g s is generated in the diagrams in two ways: 
1) due to explicit quark-gluon interaction (vertices of a hard gluon line in the diagrams 
in Figs. 10 and 11, or vertices of external gluon in the diagrams in Figs. 8 and 9), and it is 
reasonable to take it at the renormalization point fi 2 = M 2 ; 2) due to the equation of motion, 
and its normalization point should be taken in such a way that the quantity a s (0 \ ip-ip \ 0) 2 
is a renormalization group invariant. The notations in (EST) reflect this fact. 



Sum rules for Up(x) are fullfilled in wide region of x: 0.1 < x < 0.85. Borel mass M 
dependence of xu^{x) at various x is plotted in Fig. 12. One can see that it is weak in the 
whole range of x, except x < 0.1 and x > 0.75. As discussed in the Introduction the reason 
of a more strong M 2 -dependence at small and large x is connected with the fact, that our 
approach is invalid at small x and x, close to 1. It is manifested by the blow up of dimension 
6 correction at x — > and x — > 1 in eq.(^). So, the applicability domain of the sum rule can 
be found from the sum rule itself. Fig. 13 presents xu^{x) as a function of x. M 2 = 1 GeV 2 
and so = 1.5 GeV 2 , Qq = 4 GeV 2 were chosen, the other set of parameters - Aqq D and a s a 2 
is the same as in the calculation of xu w (x). 

Valence w-quark distribution in transversally polarized p-meson is given by 



xu T p (x) 



57T 



M 4 



2^P 



-X 



<p (x)E 1 : 



so 
M 2 



1 + —ln\ 
3vr ' 



01 
M 2 



K(M 2 



'(Ax 



(po{x) 
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+4/n(l - s) 



2(1 - 2x + 4s 2 )/ns" 



(f {X) 



vr 2 (0\(a s /n)G 2 \0) 
6 M 4 s 2 



2 8 -3 5 M 6 s 3 (l-s) 3N ; " 
a s (M 2 )(a s a 2 



{0\g 3 f abc GlGl x G^\0)ax) 



2 5 -3% 2 M¥(l-i) 



£(x) = -1639 + 8039a; - 15233s 2 + 10055s 3 - 624s 4 - 974s 5 



(66) 
(67) 



X(x) = 8513 - 41692s + 64589s 2 - 60154s 3 + 99948s 4 
-112516s 5 + 45792s 6 + (-180 - 8604s + 53532s 2 
- 75492s 3 - 28872s 4 + 109296s 5 - 55440s 6 )/n2 (61 



The standard value fl6j of gluonic condensate (0 | (a s /n)G 2 | 0) = 0.012 GeV 4 was taken in 



numerical calculations. Unlike the cases of u n (x) and «p(s), dimension-6 power correction, 
proportional to (0\g ' f abc G a ^ u G b v> G\ | 0), is not cancelled here. Its contribution is calculated 
with the help of the instanton gas model - eq.(^). The effective instanton radius p c was 
chosen as p c = 0.5/m. This value is between the estimations of [^(J (p c = Ifm) and 



(p c = 0.33/m). (In the recent paper |^3| the arguments were presented that the liquid gas 
instanton model overestimates higher order gluonic condensates and in order to correct this 
effect larger values of p c comparing with ETI should be used). Borel mass dependence of 



xut{x) is shown in Fig. 14. As is seen from Fig. 14, in the interval 0.2 < s < 0.65 the in- 
dependence is weak at 0.8 < M 2 < 1.2 GeV 2 . Fig.15 shows xut{x) at M 2 = 1 GeV 2 and 
Qq = 4 GeV 2 . Dashed lines demonstrate the influence of the variation of p c in the final result: 
the lower line corresponds to p c = 0.6/m and the upper - to p c = 0.4/m. Our results are 
reliable at 0.2 < s < 0.65, where d = 4 and d = 6 (separately) power corrections comprise 
less than 30% of the bare loop contribution. ((0|(a s /vr)G 2 |0) and (0|^ 3 / afec G^G^ A G^|0) 
contributions are of the opposite sign and compensate one another, a s (M 2 )a s a 2 contribution 
is negligible.) At p c = 0.4/m the applicability domain shrinks to 0.25 < s < 0.6. The 
moments of quark distributions in longitudinal p-meson are calculated in the same way, as it 
was done in the case of pion: by matching with Regge behaviour u(x) ~ l/\/s at low s and 
with quark counting rule u(x) ~ (1 — s) 2 at large s. The matching points were chosen as 
s = 0.10 at low s and s = 0.80 at large s. The numerical values of moments for longitudinally 
polarized p are 



M[ = J dxu L p (x) = 1.06 (1.05) 



M% = J xdxu L p (x) = 0.39 (0.37) (69) 
o 

The values of momenta, obtained by assuming that u(x) ~ (1 — s) at large s are given in 
the parenthesis. 
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Reliable calculation of moments in the case of transversally polarized p-meson is impos- 
sible, because of a narrow applicability domain in x and the form of w-quark distribution 
- Fig. 15, which does not allow soft matching with expected behaviour xu^(x) at small and 
large x. 

The comparison of w-quark distributions in longitudinally and transversally polarized p- 
mesons shows a strong difference of them: the curvatures have the opposite sign: u p (x) has 
a maximum at intermediate x ~ 0.5 while u^(x) has a minimum there. Strongly different 
are also the second moments in pion and in longitudinal p-meson: the total part of the 
momentum, carried by valence quarks and antiquarks — (u + d) in longitudinal p-meson is 
about 0.8, while in pion it is much less - about 0.4-0.5. The same amount as in pion one 
may expect in transverse p (Fig. 15), nevertheless, that it is impossible to calculate a precise 
number. Therefore, pion and transverse p-mesons in this aspect behave like a nucleon, where 
about 50% of total momentum is carried by gluons and sea quarks. In longitudinal p the 
situation is different - only about 20% of momentum is left for gluons and sea quarks. It 
must be mentioned, that in case of transverse p the accuracy of our results are worse than 
for longitudinal p, because the contribution of higher order terms of OPE is larger and the 
applicability domain in x is narrower. This fact, however, does not change the qualitative 
conclusion formulated above. Now let us discuss the nonpolarizied p-meson case. Quark 
distribution function u(x) in this case is equal to: 

u p (x) = (u L p (x) + 2u T p (x))/3 

and we can determine u(x) only in the region, where sum rules for u p (x) and xlF p are fulfiled, 
i.e. 0.2 < x < 0.65. In this region u p (x) is found to be very close to u^{x) (the difference in 
whole range of x is not more than 10-15%). 



5 Summary and discussion 

Let us first discuss the accuracy of our results. In case of w-quark distribution in pion 
the main uncertainty comes from the magnitude of a s (0 \ ipip | 0) 2 . For renorminvariant 
quantity (27r) 4 a s (| ipip | 0) 2 in our calculations we took the value 0.13 Gev 6 . In fact, however, 
it is uncertain by a factor of 2. (Recent determination |24[] of this quantity from r-decay 
data indicates that it may be two times larger). Also the perturbative corrections introduce 
some uncertainties, especially at large x, (x > 0.6) where the accounted LO correction is 
large. (E.g. instead of Aq C d = 200 MeV the value Aq C d = 250 MeV could be taken). 
The estimation of both effects shows, that they may result in 10-15% variation of xu n (x) - 
increasing at x < 0.3 and decreasing at x > 0.3. One should note, that our estimation of 



the second moment of quark distribution in 7r-meson (32) differ from those obtained in [pq] . 
The origin of this discrepancy is not completely clear now: partly it could be related with 
the difference of the value of parameters used in p5j (Aq^ d = lOOMeV, Ql = lGeV 2 ) from 
our choise of parameters, and partly, maybe, with the uncertainty of our estimations of the 
moments. 

In case of w-quark distribution in longitudinally polarized p-meson the uncertainties in 
a s {0 | ipip | 0) 2 do not play any role, because of higher values of M 2 and the main source of 
them is the perturbative corrections. They influence only high x domain, x > 0.5. 
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The accuracy of our results for w-quark distribution in transversally polarized p-meson 
is worse, because of a large role of d = 4 and d = 6 gluonic condensate contributions. The 
variation of xu^(x) arising from uncertainties of d = 6 gluonic condensates was shown in 
Fig.15. The gluonic condensate (0 | (a s /ir)G 2 | 0) is also uncertain by a factor 1.5. It may 
result in 30-40% variation of xwj at x ~ 0.3 — 0.4, but much less at x ~ 0.5 — 0.6. At least, 
these uncertainties do not influence the shape of w-quark distribution. (The LO perturbative 
corrections are no more than 20% at small x and negligible at large x). 

Fig. 16 gives the comparison of valence w-quark distributions in pion, longitudinally and 
transversally polarized p- mesons. The shapes of curves are quite different, especially of 
xu^(x) in comparison with xu L p [x) and xu^x). Any uncertainties in our calculations cannot 
influence this basic conclusion. The values of moments of quark distributions are also different 
- see the discussion at the end of Sec. 4. 

The main physical conclusion is: the quark distributions in pion and p-meson have not to 
much in common. The specific properties of pion, as a Goldstone boson manifest themselves 
in different quark distributions in comparison with p. SU (6) symmetry, probably, may take 
place for static properties of 7r and p, but not for their internal structure. This fact is not 
surprising. Quark distributions have sense in fastly moving hadrons. However, SU(Q) sym- 
metry cannot be selfconsistently generalized to relativistic case [26[]. We have no explanation, 
why it-quark distributions in pion and nonpolarized p-meson at 0.2 < x < 0.65 are close to 
one another - is it a pure accident or there are some deep reasons for it. 
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Figure 1: Integration region in the double dispersion representation. 
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Figure 2: Example of the non-diagonal transition 
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Figure 3: Diagram, corresponding to the unit operator contribution. Dashed lines with 
arrows correspond to the photon, thick solid - to hadron current 
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Figure 4: Examples of the non-loop diagrams of dimension 4. Wave lines correspond to 
gluons, dot means derivative, other notations as in Fig. 3. 
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Figure 5: Quark distribution function in pion, noted "1". For comparison the fit from [6], 
labelled " GRV" , is shown. 
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Figure 6: Diagrams, corresponding to the d — 4 operator contribution. Dashed lines with 
arrows correspond to the photon, thick solid - to the hadron current, wave lines correspond 
to the external gluon field 
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Figure 7: Dimension 6 diagrams corresponding to evolution of non-loop diagrams. All 
notations as in Fig. 4 



27 




+ all ombination 



Figure 8: Diagrams of dimension 6, see text. All notations as in Fig.4 
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Figure 9: Diagrams of dimension 6. External lines with dots correspond to derivatives 
external fields. All notations as in Fig. 6 
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Figure 10: Diagrams of dimension 6, corresponding to the quark propagator expansion. All 
notations as in Fig. 4 



30 



A-\ A f\ A ^ 



+ all combinations 



Figure 11: Diagrams of dimension 6, corresponding to the quark and gluon propagator 
expansion. All notations as in Fig. 4 
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Figure 12: Borel mass dependence of the quark distribution function xu^{x) at various x 
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Figure 13: quark distribution function xu 
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Figure 14: Borel mass dependence of the quark distribution function xu^(x) at various x 
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Figure 15: xu T (x) for transversally polarizied p-meson at three choices of instanton radius 
p = 2, 2.5, 3 GeV^ 1 (correspondingly, curves are labelled by r = 2, r = 2.5, r = 3) 
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Figure 16: xu(x) for p T -(curve is labelled by ro-tr), p L (is labelled by ro-1) and 7r-meson (pi) 
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